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1.  INTRODUCTION 


In  X-ray  powder  diffraction  it  is  frequently  necessary  to  predict  the 
positions  of  the  Bragg  reflections  which  correspond  to  a  given 
crystallographic  unit  cell.  This  is  a  trivial  problem  when  the  position  of 
only  a  single  Bragg  reflection  of  known  Miller  indices  is  considered;  while 
using  a  computer  it  is  relatively  simple  to  compute  the  position  of  all  the 
possible  Bragg  reflections  which  correspond  to  a  given  unit  cell  (provided  the 
computer  store  is  big  enough).  However,  it  is  not  so  easy  to  restrict  the 
calculation  to  reflections  which  are  crystallographically  distinct,  i.e. 
making  due  allowance  for  the  symmetry  properties  of  the  Laue  group  to  which 
the  unit  cell  belongs,  and  moreover  to  exclude  those  reflections  which  are  not 
allowed  by  the  general  space  group  symmetry  of  the  unit  cell  considered. 

FIREBIRD3  has  been  written  to  do  all  these  things  and  to  execute  calculations 
which  relate  to  crystal  densities  and  interplanar  angles.  Moreover,  the 
program  can  vary  a  given  set  of  lattice  parameters  over  a  pre-determined  range 
and  since  it  can  also  compare  the  calculated  positions  of  Bragg  reflections 
with  those  actually  observed  it  can  be  used  as  a  means  to  index  X-ray  powder 
diffraction  patterns.  The  program  can  also  produce  blocks  of  Miller , indices 
in  a  form  suitable  for  transcription  as  input  to  the  program  FIREFLY^  '  and  in 
this  sense  should  be  considered  as  complementary  to  it. 

FIREBIRD  represents  one  of  a  series  of  programs  developed  at  the  UKAEA  .  . 

Northern  Research  Laboratories  (Springf ields)  (SL)  and  at  UXAEA, ,Risley . 

Various  programs  have  been  developed  elsewhere  along  these  lines ^  '  but 

FIREBIRD3  is  believed  to  represent  the  most  satisfactory  program  at  the 
present  time. 

FIREBIRD3  is  one  of  a  suite  of  programs  developed  at. SL_.for  the  solution  and 
interpretation  of  X-ray  powder  diffraction  records,  ’  ’  ’  '  and  as  such  it 

also  carries  the  designation  HELLFIRE7.  A  complete  list  of  the  programs  in 
the  HELLFIRE  suite  is  given  in  Appendix  8. 

2.  EQUATIONS  SOLVED 

Only  standard  crystallographic  equations  have  been  used.  They  are  summarised 
here,  but  given  in  full  in  Appendix  1.  The  unit  cell  volume,  V,  is  calculated 
from  the  lattice  parameters  (a,  b,  c,  a,  g  and  y)  of  the  crystallographic  unit 
cell  using  the  general  formula  for  a  triclinic  unit  cell.  '  The  reciprocal 
lattice  parameters  are  derived  from  the  lattice  parameters  given  by  standard 
formulae.  The  number  of  lattice  points,  n,  within  the  limiting  (Ewald) 

sphere  is  given  by: 

n  =  32trV/3A3 

where  A  is  the  first  X-ray  wave-length  considered.  The  'FLINT  CODES1^^ 

which  will  be  described  later,  is  a  set  of  numbers  which  specifies  the  general 
conditions  limiting  the  possible  reflections  for  a  particular  space  group. 

The  Bragg  angle, 8*. has  been  calculated  using  the  formula  appropriate  to  each 
crystal  system.  '  The  interplanar  spacing,  d(X),  is  given  from  the  Bragg 
equation  □ 

□ 

A  =  2d  sin  9  - - 


where  the  order  of  the  reflection  has  been  ignored,  since  it  is  conventional  _ 

to  incorporate  it  within  the  Miller  indices  h,  k  and  i  which  describe  a 
reflecting  plane.  The  program  never  refers  to  the  hexagonal  Miller  index  i,  cJ(le 
which  equals  -  (h+k),  while  it  refers  to  the  rhorobohedral  Miller  indices  as 

id/or 

cist  i  jpocial 
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h,  k  and  l  rather  than  p,  q  and  r.  The  program  uses  'NSYS  codes'  to  define 
the  various  Laue-symmetry  groups.  .These  codes,  given  in  Appendix  2,  were 
developed  for  the  program  FIREFLY.  '  They  are  used  to  derive  the 
multiplicity  of  each  reflection  output  and  to  define  the  non-equivalent 
arrangements  of  each  Miller  index.  Interplanar  angles  relative  to  one  or  more 
given  crystallographic  reference  planes  are  calculated  using  the  general 
formulae  for  a  triclinic  unit  cell.  1  The  Nelson-Riley  function  referred  to 
in  certain  outputs  is  given  by  the  expression: 


l  f  cos  2  8  cos  2  9^1 

t  {  sin  8  8  J  , 

( 19) 

it  is  used  in  certain  crystallographic  extrapolation  procedures. 

The  program  calculates  the  theoretical  ('crystal')  density,  p,  from  the 
volume,  V,  of  the  unit  cell,  the  molecular  weight  of  its  contents,  M,  and  the 
Avogadro  number.  A,  using  the  relation 

p  =  M/AV 


3.  GENERAL  AND  SPECIFIC  FEATURES 

FIREBIRD3  uses  the  standard  methods  of  data  input  applicable  to  FORTRAN  77; 
the  detailed  form  of  the  data  input  is  described  later.  The  program  operates 
on  the  equations  referred  to  above.  Thus,  for  a  given  set  of  unit  cell 
dimensions,  FIREBIRD  calculates  and  outputs  the  unit  cell  volume,  reciprocal 
lattice  parameters,  and  the  number  of  reciprocal  lattice  points  corresponding 
to  the  first  wavelength  considered  within  the  limiting  (Ewald)  sphere, 
(although  the  program  makes  no  further  use  of  the  latter  information).  Then, 
provided  it  is  given  data  which  enable  it  to  determine  the  molecular  weight  of 
the  contents  of  the  unit  cell,  it  infers  the  crystallographic  density. 

The  main  role  of  the  program  is,  however,  to  compute  a  list  of  the  Bragg 
reflections,  subject  to  certain  constraints,  for  all  possible  Miller  indices 
which  can  be  produced  by  a  given  crystallographic  unit  cell,  using  X-ray  (or 
electron)  powder  diffraction  techniques.  It  does  this  by  calculating  sin28 
from  the  equation  appropriate  to  the  crystal  system  and/or  setting,  and/or 
lattice  of  the  unit  cell  and  the  appropriate  wavelength  being 
considered.  * 

Various  X-ray  wavelengths,  Cu  Ka.1,  Cu  Ka2,  Cu  Ka  (weighted  mean), 

Cu  K(5,  Cr  Kaj ,  Cr  Ka2,  Cr  Ka  (weighted  mean),  Fe  Ka2 ,  Fe  Ka2 , 

Fe  Ka  (weighted  mean),  Co  Ka2 ,  Co  Ka2,  Co  Ka  (weighted  mean),  Mo  Kalf  Mo  Ka2 , 
Ag  Ka2 ,  and  Ag  Ka2  are  stored  within  the  program  and  referred  to  by  codes  in 
the  data  input.  All  this  information,  which  is  detailed. in  Appendix  3,  has 
been  taken  or  derived  from  Henry,  Lipson  and  Wooster.  ’  Since  these 
wavelength  values  are  subject  to  revision  in  the  literature  the  preset  values 
can  be  overwritten,  as  required,  in  the  data  input. 

Provision  is  made  for  the  program  to  make  systematic  variations  in  the  lattice 
parameters  given,  so  altering  the  predicted  reflection  positions.  When  these 
data  are  output  alongside  data  which  refer  to  observed  reflection  positions, 
it  becomes  possible,  in  principle,  to  assign  Miller  indices  to  any  given  set 
of  observations.  However,  as  the  symmetry  of  the  unit  cell  decreases,  the 
computer  time  required  rapidly  becomes  unreasonable.  In  the  program's 
calculations,  the  crystal  system  is  either  derived  directly  from  the  Laue 
group  input  in  coded  form,  or  inferred  from  the  space  group. 
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The  next  section  describes  in  some  detail  how  the  Miller  indices  are  derived 
for  each  crystal  system;  while  a  further  section  explains  how  the  general 
conditions  which  govern  the  allowed  reflections  from  a  particular  space  group 
are  operated. 

The  number  of  reflections  listed  is  limited  to  a  maximum  of  1000,  or  some 
other  lower  predetermined  limit,  or  up  to  a  preset  Bragg  angle.  The 
reflections  are  output  in  order  of  increasing  angle  irrespective  of  the 
wavelength,  which  is  referred  to  in  each  case  by  a  numerical  code. 

Two  statements  appear  at  the  end  of  each  example.  The  first:  THE  NEXT  HLK  TO 
BE  CONSIDERED  WAS  ...  indicates  the  nex  highest  value  of  either  h,  k  or  l 
beyond  that  considered  in  the  output  as  .'uld  have  been  expected  according  to 
the  logic  applicable  to  the  crystal  system  being  considered.  The  second 
statement  THERE  WERE  . . .  HKL  COMBINATIONS  AVAILABLE  ...  indicates  the  number 
of  valid  sets  of  Miller  indices  for  a  particular  example  between  0  and  90°9 
irrespective  of  those  actually  output. 

3.1  THE  CONDITIONS  WHICH  DETERMINE  THE  GENERATION  OF  h,  k  AND  l  VALUES  FOR 
EACH  LAUE  GROUP 

The  crystal  system  is  either  derived  from  the  Laue  group,  input  in  coded  form, 
or  inferred  from  the  space  group,  subject  to  the  fact  that  the  program  is  so 
written  that  the  Laue  group  given  overrides  the  space  group.  It  will, 
moreover,  be  seen  from  the  input  specification,  that  various  details  of  the 
space  group  or  Laue  group  can  be  suppresssed  in  the  data  output. 

In  the  cubic  system,  the  maxium  possible  value  of  the  Miller  index,  h,  is 
first  inferred,  on  the  basis  that  sin29  cannot  exceed  unity,  or  a  value 
coresponding  to  the  preset  value  of  9;  it  is  a  value  which  k  and  l  cannot 
exceed.  Every  wavelength  is  considered  in  this  calculation.  All  possible 
positive  arrangements  of  h,  k  and  l  are  then  computed  subject  to  the 
conditions  that  for  the  Laue  group  m3m,  h  ;  k  M;  while  for  the  Laue  group 

m3,  the  m3m  list  is  further  extended  by  the  interchange  of  h  and  k,  i.e.  k  can 

exceed  h. 

In  the  tetragonal  crystal  system  maximum  values  of  h  and  i  are  deduced  as 
above.  For  the  Laue  group  4/mmm,  all  possible  positive  non-identical 
arrangements  of  h,  k  and  t  are  next  inferred,  subject  to  the  overriding 
condition  that  k  cannot  exceed  h.  The  tetragonal  Laue  group  A/m  is  treated  as 

for  A/mmm  except  that  k  can  become  negative,  unless  it  equals  h. 

For  the  orthorhombic  system,  Laue  group  mmm,  each  maximum  value  of  h,  k  and  i 
is  first  inferred.  The  program  then  evaluates  all  different  possible  positive 
arrangements  of  h,  k  and  !.. 

For  the  hexagonal  crystal  system  maximum  values  of  h  and  l  are  inferred  as 
above.  For  the  Laue  group  6/mmm  all  possible  positive  non-identical 
arrangements  of  h,  k  and  l  are  then  inferred  subject  to  the  overriding 
restriction  that  k  cannot  exceed  h.  However,  for  the  Laue  group  6/m  k  can 
additionally  exceed  h  except  for  hOS.  when  OhS.  is  not  considered.  When  the 
trigonal  Laue  group  3  is  considered  as  hexagonal  it  is  treated  as  for  Laue 
group  6/m  above  except  that  l  can  also  take  negative  values,  i.e.  hki,  kh4, 
hk£  and  khU  are  all  considered.  Similarly  when  the  Laue  grou£  3m  is  treated 
as  hexagonal  it  is  evaluated  as  for  3  except  that  for  hhl,  hhl  is  not 
considered  and  for  hkO,  khO  is  not  considered. 
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For  the  monoclinic  crystal  system,  Laue  group  2/m,  second  setting  (the 
International  Union  of  Crystallography  preferred  setting),  the  maximum  values 
of  h,  k  and  2.  are  computed;  the  mixed  index  term  (see  Appendix  1)  is  ignored 
and  each  maximum  value  is  arbitrarily  tripled  to  allow  for  this.  All  possible 
positive  non-identical  arrangements  of  k  and  2.  are  then  considered  with  all 
possible  positive  and  negative  values  of  h,  except  that  when  2.  is  zero  h 
cannot  be  negative.  The  first  setting  of  the  monoclinic  crystal  system  is 
treated  similarly  but  k,  instead  of  h,  can  become  negative,  except  when  h  is 
zero.  (If  for  the  first  setting  a  =  a,  b  =  b,  c  =  c,  y  =  y,  then  if  the 
second  setting  is  applied  to  the  same  crystal  structure  a'  =  a,  b'  =  c,  c'=  b, 
(3  »  ( 180°  -  y). 

In  the  trigonal  crystal  system  considered  with  a  rhombohedral  lattice  and  for 
the  3m  Laue  group,  nested  loops  are  applied  wherein  h  is  increased,  followed 
by  k  and  then  2,  such  that  these  values  take  only  positive  values  and  2.  does 
not  exceed  k,  and  k  does  not  exceed  h.  All  possible  values  of  h,  k  and  l  are 
considered  until  either  1000  arrangements  have  been  found  which,  allowing  for 
subsequent  tests,  can  be  output,  or  1000  arrangements  have  been  considered,  or 
the  data  input  have  imposed  some  other  limit.  In  this  way  the  (h2  +  k2  +  2.2) 
term  is  found  for  the  standard  rhombohedral  formula,  see  Appendix  1.  The 
(hk  +■  kH  4-  hi)  term  is  next  derived  as  follows:  the  values  of  h,  k  and  2. 
already  found  are  first  considered  as  positive,  each  hk,  ki  and  hi  term  is 
then  allowed  separately  to  become  negative  so  that,  in  all,  four  values  are 
found  for  the  (hk  +  ki  +  hi)  term.  However  only  different  values  of  this  term 
are  considered  further_for  each  set  of  values  of  h,  k  and  i,  i.e.  432  does  not 
exclude  520.  For  the  3  Laue  group  treated  as  rhombohedral  the  same 
considerations  apply  except  that  for  hki,  hik  must  also  be  output. 

For  the  triclinic  crystal  system,  Laue  group  I,  nested  loops  are  applied  in 
which  h  is  increased,  followed  by  k,  followed  by  i.  Initially  h,  k  and  t  are 
all  treated  as  positive  and  i  }  k  }  h.  Next  all  combinations  of  each  of  these 
hki  arrangement  are  considered  in  which  h,  k  and  i  are  independently  either 
positive  or  negative,  thus  there  are  72  possible  arrangements  based  on  321. 
These  computations  stop  when  ten  thousand  sets  of  Miller  indices  have  been 
considered,  unless  either  a  thousand  values  of  0  have  been  output,  or  some 
other  limit  has  been  introduced  in  the  data  input. 

3.2  THE  GENERAL  CONDITIONS  LIMITING  THE  POSSIBLE  REFLECTIONS  FROM  A  SPACE 
GROUP  -  THE  FLINT  CODES 

At  this  stage  the  program  has  calculated  the  angular  position  of  all  the  Bragg 
reflections  corresponding  to  a  particular  crystal  system  within  certain 
limits.  However,  considerable  numbers  of  reflections  may  be  forbidden  by  the 
general  symmetry  of  a  particular  space  group  e.g.  for  a-alumina  only  eighty 
reflections  are  actually  allowed  by  the  trigonal  space  group  167  whereas  in 
terms  of  the  crystal  system  alone  270  reflections  would  be  predicted.  For 
this  reason  nine  digits,  the  FLINT  codes,  have  been  derived  which  summarise 
the  reflections  allowed  by  a  particular  space  group.  These  codes  are  input  to 
the  program  as  the  integers  NF1  to  NF9,  and  are  referred  to  within  the  program 
as  NSYM1  to  NSYM9 . 

The  development  of  these  codes  is^b|sed  directly  on  Table  4.1.9  and  Section 

4.3  of  the  International  Tables.  Table  4.1.9  was,  however,  specifically 

intended  for  X-ray  diffraction  studies  of  single  crystals,  and  some 
modifications  have  been  necessary  in  the  derivation  of  the  NSYM  codes  to  make 
them  meaningful  in  X-ray  powder  diffraction  calculations.  These. modifications 
are  noted  below.  Since  the  FLINT^ggdes  were  compiled  by  Stocks'  ’  it  has 
become  necessary  to  update  them. 
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The  implications  of  the  nine  FLINT  codes  or  NSYM  codes  are  described  fully  in 
Appendix  4.  There  the  NSYM  codes  are  listed  together  with  the  conditions 
which  the  Miller  indices  h,  k  and  i  must  satisfy  for  the  reflection  to  be 
allowed.  The  relevant  portions  of  the  short  Hermann-Maugain  space  groups  are 
also  listed.  The  full  space  group  symbols  are  not  considered  as  these  lead  to 
'redundant'  conditions,  i.e.  conditions  already  catered  for  by  the  symmetry 
properties  of  the  short  symbol.  Some  features  of  these  symbols  are  not 
strictly  consistent  between  crystal  classes,  thus  in  the  orthorhombic  space 
group  Pca2i  (No.  29),  2X  refers  to  a  [001)  screw  axis,  whereas  in  the  ,  . 
tetragonal  space  group  14^2  (No.  98),  4!  refers  to  a  [001]  screw  axis. 

In  broad  outline,  the  operation  of  the  NSYM  codes  is  as  follows: 

NSYM1  covers  the  conditions  imposed  by  the  Bravais  lattices,  i.e.  face  and 
body  centring,  and  various  types  of  rhombohedral  and  hexagonal  indexing. 
Conditions  8  and  10  do  not  apply  directly  to  the  space  groups  listed  in  the 
International  Tables  but  they  are  included  to  present  a  complete  picture;  the 
program  allows  them  to  be  input  if  required. 

NSYM2  covers  the  restrictions  imposed  on  reflections  of  the  type  Okt  by  the 
presence  of  {100}  glide  planes.  Strictly  it  should  apply  not  only  to  the 
orthorhombic  systems  but  to  the  tetragonal  and  cubic  systems.  However,  by 
convention.  X-ray  powder  reflections  for  crystals  of  the  cubic  systems  are 
described  as  hkO  rather  than  Oki  and  so  this  type  of  reflection  is  considered 
under  NSYM4.  Similarly  since  FIREBIRD2  refers  to  hOH  rather  than  its 
equivalent  OkH.  in  the  tetragonal  system,  reflections  of  this  type  are 
considered  under  NSYM3  rather  than  under  NSYM2. 

NSYM3  considers  the  implications  of  (010)  glide  planes  in  a  structure  thus  it 
considers  reflections  of  the  type  hOt  in  the  monoclinic  (second  setting), 
orthorhombic  and  tetragonal  (see  NSYM2)  systems. 

NSYM4  treats  reflections  of  the  type  hkO  and  refers  to  the  monoclinic  (first 
setting),  orthorhombic,  tetragonal  and  cubic  (see  NSYM2)  crystal  systems. 

NSYM4  deals  with  the  presence  of  (001)  glide  planes  in  a  structure. 

NSYM5  covers  restrictions  imposed  on  reflections  of  the  type  hhi  and  applies 
to  the  tetragonal,  trigonal,  hexagonal  and  cubic  systems.  NSYM5  covers  the 
presence  of  {110}  glide  planes  in  a  structure.  Under  NSYM5,  code  2  covers  not 
only  the  condition  that  t  must  be  even,  but  also  the  n  glide  condition  which 
requires  that  (2h  +  l)  must  be  a  multiple  of  four. 

NSYM6  deals  with  reflections  of  the  type  hOO  and  treats  the  orthorhombic, 
tetragonal  and  cubic  systems.  It  covers  the  effects  of  [100]  screw  axes  in  a 
structure. 

NSYM7  covers  the  restrictions  imposed  on  reflections  of  the  type  OkO  by  [010] 
screw  axes  in  the  monoclinic  (second  setting)  and  orthorhombic  systems. 

NSYM8  treats  the  implications  of  [001]  screw  axes  on  reflections  of  the  type 
00J.  in  the  monoclinic  (first  setting),  orthorhombic  and  tetragonal  systems  and 
z  screw  axes  in  the  hexagonal  systems. 

NSYM9  covers  restrictions  imposed  on  reflections  of  the  type  hhl  (encountered 
in  FIREBIRD2  as  hO!,)  by  {1120}  glide  planes  in  the  trigonal  and  hexagonal 
systems . 
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The  values  of  NSYM  corresponding  to  each  space  group  are  given  in  Appendix  5, 
and  this  table  is  accessed  via  the  program  NSYM  from  a  disc  file  to  the  main 
program  FIREBIRD3. 

3.3  PROGRAM  LISTING  AND  SUBROUTINE  FUNCTIONS 

The  complete  listing  of  the  program  FIREBIRD3  is  given  on  micro-fiche  (inside 
front  cover).  The  role  of  the  various  constituent  subroutines  is  summarised 
in  Appendix  6.  The  actual  instructions  for  compiling  and  running  FIREBIRD3  on 
the  UKAEA's  Risley  ICL  3980  computer  are  given  in  Appendix  7. 

4.  FORM  OF  INPUT 

Data  are  input  to  the  program  from  a  standard  72  character  record  data  file 
(in  lieu  of  standard  80-column  IBM  cards).  In  general  only  one  variable 
appears  on  a  record  (a  single  IBM  card)  except  where  stated  below.  With  the 
exception  of  two  variables,  MAXNUM  and  NUMOBS,  which  are  free  format,  all  the 
variables  are  input  in  mixed  format.  In  the  input  lists  given  below  for 
FIREBIRD3  many  variables  are  enclosed  in  square  brackets.  Whether  or  not  such 
variables  are  input  depends  upon  the  program  options  specified.  In  the  format 
specifications  given  below,  F  refers  to  a  real  number,  I  to  an  integer  number, 
10  to  a  free  format  integer  and  DO.O  to  a  floating  point  number.  The  lists 
which  define  the  order  of  variables  input,  the  names  of  the  variables  within 
the  program,  their  respective  formats,  range  of  meaningful  values  assigned  to 
each  variable  and  their  respective  roles  are  given  below. 


Name  of  Inclusive  range  of 

input  Format  values  Role  of  variable 

variable 


NO  13 


TITLE 


0UTPT1  II 


0UTPT2  II 


ND-R-91KS) 


1  to  100 


6  records 


2 


0 

2 


This  record  appears  once  only  at 
the  beginning  of  a  set  of 
examples.  It  defines  the  number 
of  examples  to  be  considered. 

These  six  records  must  appear  for 
each  example,  although  some  or 
all  of  them  may  be  blank.  They 
contain  up  to  6  x  80  characters 
from  the  256  character  EBCDIC  set 
which  are  necessary  to  define  a 
'title'  for  each  example. 

131  column  wide  output  required 
corresponding  to  calculated  and 
observed  Bragg  reflection 
positions. 

No  output  required  which  relates 
to  Bragg  reflection  positions. 

Output  of  interplanar  angles, 
see  the  variable  BALL  which 
appears  later,  required. 

N.B.  0UTPT2  set  at  2 
overrides  all  settings  of  0UTPT1 . 
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Name  of 
input 
variable 

Format 

Inclusive  range  of 
values 

Role  of  variable 

0UTPUT2 

cont 

1 

Outputs  blocks, of  Miller  indices 
in  a  FIREFLY  compatible  form. 

0 

No  0UTPT2  option  specified 

OBTHET 

12 

3 

The  output  data  include  the 

1.  1  .  loRln  d,  sin  8,  sin  8. 
d  d2  X 

the  Nelson  Riley^^  function, 
and  cos28;  no  observations  are 
given. 

2 

Observed  Bragg  angles  given  in 
°28 . 

1 

Observed  Bragg  angles  given  in 

“9. 

0 

No  observed  data  given. 

-1 

Observed  Bragg  angles  given  as  d, 
interplanar  spacings,  in  A. 

-2 

Observed  Bragg  angles  given  as 
sin29  values. 

NSYS 

13 

15  to  2 

Code  defining  the  Laue  symmetry 
group  (see  Appendix  2),  NB  NSYS 
overrides  SG  (see  below)  in 
program. 

1 

Unallocated  value. 

0 

The  program  infers  NSYS  from  SG 
and  SETT  (see  later). 

-1 

When  NSYS  is  negative,  NF1  to  NF9 
are  read  in  (see  below)  but  when 

NSYS  is  -1,  LAUE  SYMMETRY  GROUP 
UNDEFINED  is  output.  A  'dummy' 
value  of  SG  must  be  input  later, 
so  that  the  necessary  lattice 
parameters,  'A'  etc.  (see  below) 
will  be  read  in.  In  general  SG 
would  then  be  negative. 

-2  to  -15 

NF1  to  NF9  are  read  (see 

above/below).  These 
negative  NSYS  values  are  treated  as 
their  positive  equivalents  but  no 
knowledge  is  required  of  the 
space  group  and  SG  may  be  read 
in  later  as  0  (when  it  is  not 
referred  to  in  the  output). 


ND-R-911(S) 


9 


Inclusive  range  of 
values 


Role  of  variable 


Name  of 

input  Format 
variable 


[NF1] 

[NF2] 

[NF3] 

[NF4] 

[NF5]  912 

[NF6] 

[NF7] 

[NF8] 

[NF9] 


SG  14 


SETT  II 


1  to  10 
1  to  5 
1  to  5 
1  to  5 
1  to  3 
1  to  3 
1  to  2 
1  to  5 
1  to  2 


230  to  1 
0 


-1  to  -230 


6 


3 


2 


1 


These  codes  (given  on  a  single 
record)  are  only  read  when 
NSYS  is  negative.  They  define  the 
conditions  which  limit  the  possible 
reflections  from  crystal  system,  see 
Appendix  2.  Except  when  NSYS  is 
negative  these  codes  are  referred  to 
within  the  program  which  has  access 
to  NF1  to  values  of  NF9  matched 
against  each  spacegroup  (SG). 

The  space  group  applicable. 

No  conditions  are  applied  as  to 
whether  or  not  a  reflection  is 
allowed.  If  SG  is  0,  NSYS 
cannot  be  zero. 

As  for  SG  equals  1  to  230  but 
Space  group  is  not  referred  to 
in  the  output. 

NB.  NSYS  always  overrides  SG 
and  SETT. 

Defines  a  crystal  setting  or 
axes  thus: 

Refers  the  trigonal  space  groups 
143  to  167  to  hexagonal  axes;  it 
should  be  used  for  all  hexagonal 
examples . 

Refers  the  trigonal  space  groups 
143  to  167  to  rhombohedral  axes, 
where  this  is  valid. 

Refers  the  monoclinic  space 
groups  3  to  15  to  their  second 
setting  (y-axis  unique;  BETA  - 
see  later  -  required  in  input) 

(the  preferred  I.U.Cr  setting). 

Refers  the  monoclinic  space 
groups  3  to  15  to  their  first 
setting  (z-axis  unique;  GAMMA  - 
see  later  -  required  later) 

(this  setting  is  not  recommended 
by  the  I.U.Cr). 
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Name  of 

Inclusive  range  of 

input 

variable 

Format 

values 

Role  of  variable 

SETT 

0 

Space  groups  3  to  15  are  referred 

Cont. 

to  their  second  setting  (see 
above).  Space  groups  143  to  194 
are  referred  to  hexagonal  axes 
(see  above).  For  all  other  space 
groups  SETT  is  read  in,  but 
ignored  by  the  program. 

REV 

11 

Defines  the  setting  used  in  the 
trigonal  space  group  when 
hexagonal  axes  are  used: 

0 

Obverse  setting  applies. 

1 

Reverse  setting  applies. 

A 

The  crystal  lattice  parameters 

[B] 

3F8.4 

999.9999 

a,  b,  c,  a,  8,  y  are  given  on 

[C] 

a  single  record;  the 
particular  ones  input  depend  upon 
the  cyrstal  system  e.g.  the 

[ALFA] 

180.0000 

trigonal  (rhombohedral  setting) 

[BETA] 

3F9.4 

-180.0000 

requires  A  and  ALFA.  The  angles 

[GAMMA] 

are  given  in  degrees. 

NUMWAV 

12 

13 

The  number  of  wavelengths  used 
(the  results  are  sorted  in  order 
of  increasing  Bragg  angle). 

WAVECO 

13 

1  to  17 

List  of  wavelength  codes 
applied,  see  Appendix  3. 

0 

Unallocated  value  in  program. 

-1  to  -17 

In  the  example  given  the  program 
values  of  the  X-ray  wavelengths 
corresponding  tc  all  the  positive 
WAVECO  values  are  replaced  by  the 
following  WAVECO  values  of 

WAV ( N ) . 

(WAVN(J)] 

F11.8 

99.99999999 

List  of  wavelengths  to  be  input 

(see  WAVECO  above). 
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Name  of 

input  Format 
variable 


Inclusive  range  of 
values 


Role  of  variable 


CUTOFF  F9 . 5 


1  x  10' 5  to  180.00000 


1  x  10‘ 5  to  90.00000 


999.99999  to  1  x  10‘ 5 


1  x  10- 5  to  1.00000 


This  value,  which  overrides  MAXNO 
(see  below)  unless  set  at  zero 
defines  a  limiting  value  beyond 
which  no  data  are  output,  although 
the  program  actually  manipulates 
values  beyond  this  limit. 

corresponding  to  OBTHET  set  at  2 
(observations  in  °29  input). 

corresponding  to  OBTHET  set  at  1 
or  0  (observations  in  ’8  input  or 
no  observations  input). 

corresponding  to  OBTHET  set  at  -1 
(observations  input  as  d 
spacings) . 

corresponding  to  OBTHET  set  at  -2 
(observations  input  as  sin 2 6 
values ) . 


0  to  1000 


MAXNUM  10  any  positive  integer 


NUMOBS 


S  1000 


The  maximum  number  of  Bragg 
reflections  listed  in  order  of 
ascending  Bragg  angle  in  the 
output.  This  value  is  not 
allowed  to  exceed  1000,  although 
up  to  10000  sets  of  Miller 
indices  may  be  considered  by  the 
program.  NB.  MAXNO  is  over¬ 
ridden  by  CUTOFF  so,  for  example, 
if  MAXNO  is  0  it  is  CUTOFF  that 
decides  the  size  of  the  output. 

This  value  refers  to  the  maximum 
number  of  increments  allowed  when 
the  lattice  parameters  are  varied 
incrementally  (see  main  text). 

If  MAXNUM  is  0  it  defaults  to  20. 
NB.  An  incautious  use  of  MAXNUM 
could  prove  very  expensive! 

Data  which  compare  observed  and 
calculated  information  are  only 
output  when  the  number  of  "atched 
data  pairs  within  TOL  (see  Delow) 
is  not  less  than  NUMOBS . 


Observed  and  calculated  data  are 
considered  to  be  matched  when 
they  differ  by  less  than  TOL  in 
sin29 .  (If  OBTHET  is  0  then  TOL 
is,  of  course,  0). 
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Name  of 
input 
variable 

Format 

Inclusive  range  of 
values 

Role  of  variable 

RHO 

DO .  0 

0 

Continue  reading  data  at  OBTH(J), 
see  below. 

1 

Uses  Avogadro's  number  set  within 
program  at  6.022521  x  1023. 

1.0 

Resets  Avogadro's  number  to  the 
value  assigned  to  RHO. 

[TOAT] 

12 

99 

Number  of  types  of  atom  within 
the  unit  cell. 

[ATWT(I) ]] 
[NOAT(I) ] J 

F9.5 

F8.3 

(999.99999 

[9999.999 

Atomic  weight  of  Ith  atomic 
species  (on  the  same  record  as) 
the  number  of  the  Ith  type  of 
atom  present  within  the  unit  cell 
(note  this  can  be  fractional  to 
allow  for  different  isotopic 

NB.  I  equals  TOAT  (see  above). 

OBTH( J ) 1 
BWAV ( J ) j 

F9.5 

13 

f999. 99999 
[17  to  -17 

List  of  matched  pairs  (on  single 
records)  of  observed  Bragg 
reflection  positions  (OBTHET 

above)  and  a  wavelength  code 
corresponding  to  WAVECO  (see 
above ) . 

NB.  This  list  must  be  in  order 
or  ascending  9.  When  BWAV(J)  is 
set  at  zero,  it  and  all 
subsequent  observations  are 
referred  to  the  last  wavelength 
code  read  in,  and  BWAV(J)  need 
not  be  input  again. 

OBTH(J)  BWAV(J)  set  at  zero 
terminates  the  reading  of 
[OBTH(J)]  [BWAV(J)]  and  data 
input  recontinues  at  DELA. 

[BALL]  13  999  If  0UTPT2  is  set  at  2,  BALL  is 

read  here.  It  corresponds  to 
the  number  of  Miller  indices 
following:  however  the  program 
will  only  compute  data  which 
refer  to  a  maximum  of  ten  sets  of 
reference  MILLER  INDICES. 

[HR(IB)]  314  999  The  number  of  records  following  is 

[ KJR( IB ) ]  given  by  BALL  (above).  These 

(LR(IB)]  records  contain  values  of  the 

Miller  indices  of  the  reference 
planes  (see  main  text).  The 
program  reads  all  these  records 
but  only  considers  the  first  ten 
of  them. 
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Name  of 

input  Format 
variable 


Inclusive  range  of 
values 


Role  of  variable 


DELA 

F6.4 

[DELB] 

5F6.4 

[DELC] 

or  as 

[DELALF] 

appro¬ 

[DELBET] 

priate 

[DELGAM] 

[MAXA]  6F8.4 

[MAXB]  or  as 

[MAXC]  appro- 

[MAXALF]  priate 

[MAXBET] 

[MAXGAM] 


999.9999  to 
001.0000 
-180.0000  to 
0180.0000 


4.1  RESTRICTIONS  ON  INPUT 


If  DELA  is  1.0000,  the  program 
moves  on  to  the  TITLE  of  the  next 
example  otherwise  it  also  needs 
DELB  etc.,  all  on  another  record. 
DELA,  DELB  etc.,  correspond  to 
the  respective  increments  which 
will  be  applied  to  the  variables 
A,  B,  C,  ALFA,  BETA  and  GAMMA, 
but  note  that  values  of  DELA  etc. 
are  only  input  as  are  logically 
consistent  with  the  values  of 
NSYS,  SG  and  SETT  already  input. 

These  data,  all  on  one  record, 
define  the  limits  beyond  which 
the  operation  of  DELA,  DELB, 

DELC,  DELALF,  DELBET  and  DELGAM 
cannot  increase  A,  B,  C,  ALFA, 
BETA  and  GAMMA,  as  are  logically 
consistent  with  the  values  of 
NSYS,  SG  and  SETT  already  input. 


Apart  from  the  field  widths  described  above  the  only  restriction  on 
the  program  is  that  the  number  of  observations  input  cannot  exceed  1000. 

5.  DESCRIPTIONS  OF  OUTPUT 


The  number  of  reflections  output  is  limited  to  a  maximum  of  1000 
irrespective  of  wavelength  unless  the  number  of  reflections  over  the 
angular  range  considered  is  less  than  that  value.  The  maximum  angle 
considered  is  normally  90!>9  but  it  can  be  restricted  by  the  data  input. 
Subject  to  the  latter  overriding  conditions,  a  value  can  be  inserted  in  the 
data  input  which  limits  the  actual  number  of  reflections  in  rhombohedral 
and  triclinic  examples  before  truncating  the  data  for  output. 

Data  may  be  output  in  several  ways.  Each  reflection  is  always 

referred  to  in  terms  of  its  Miller  indices,  wavelength  code  number,  9,  20, 

sin29,  d  and  multiplicity.  If  required 

(a)  ,  log  d,  sin  9,  (sin  9) A  ,  the  Nelson  Riley  function  and 

d  d2  10 

cos29  may  be  output  or 

(b)  it  lists  data  which  refer  to  actual  observed  reflection 
positions  as  follows:  9,  29,  sin29,  d,  as  well  as  (calculated 
minus  observed  values)  of  8,  28,  sin28  and  d.  Each  observation 
is  matched  with  the  first  calculated  reflection  which  lies 
within  a  pre-assigned  sin28  window  (TOL). 

(c)  it  lists  a  reference  plane  and  the  interplanar  angle  between  it 
and  each  calculated  plane. 
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In  addition  blocks  of  h,  k  and  t  can  also  be  output  in  a  form  directly 
suitable  for  transcription  as  input  to  FIREFLY. 

Examples  of  the  various  forms  of  output,  corresponding  to  various 
inputs,  are  given  on  micro-fiche  as  well  as  the  FLINT  codes. 


6.  CONCLUSIONS 

The  program  has  been  tested  thoroughly  against  a  test  set  of  examples,  see 
the  micro-fiche  (inside  front  cover). 

The  program  listing  is  given  on  micro-fiche  (also  inside  front  cover),  and  the 
role  of  the  various  subroutines  is  summarised  in  Appendix  6,  Appendix  7 
lists  the  instructions  for  compiling  and  running  FIREBIRDS  on  the  Risley 
ICL  3980  computer.  This  program  represents  a  final  working  version  and  no 
significant  modifications  are  anticipated;  although  it  is  now  realized  that  in 
the  triclinic  case  the  two  equivalent  reflections  hki  and  hk£;  hk£  and  hkt, 
etc  are  all  output  by  FIREBIRD  3. 
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8.  APPENDIX  1 

Crystallographic  equations  used 

8.1  Volume  of  unit  cell,  V 

V  =  2  abc  /{sin  s.  sin(s-a).  sin(s-g).  sin(s-y)} 
where  s=(a+g+y)/2 

8.2  Volume  of  reciprocal  unit  cell,  V* 

V*  ■  a*b*c*  sina*sing*siny* 

8.3  Reciprocal  lattice  parameters  for  the  different  crystal  systems 

(i)  cubic  (NSYS  =  8  or  15) 

a*  =  b*  =  c*  =  -;  a*  =  g*  =  y*  =  90° 

(ii)  rhombohedral  (NSYS  7  or  14) 

a*  =  b*  =  c*  =  (a2  sina)/V,  where  v  =  a3  [1-3  cos2a  +  2  cos3a]^ 
cos  a*  =  cosg*  =  cosy*  ”  '  cosa/ ( 1  +  cosa) 

(iii)  hexagonal  (NSYS  =  6,  11,  12  or  13) 

a*  -  b*  =  2/a/3;  c*  -  A;  a*  =  g*  *  90°;  y*  ”  60° 

(iv)  tetragonal  (NSYS  =  5  or  10) 

a*  =  b*  =  ~  c*  =  a*  =  g*  =  y*  =  90° 

(v)  orthorhombic  (NSYS  =  4) 

a*  =  I;  b*  =  i;  c*  =  — ;  a*  =  g*  =  y*  =  90° 
a  d  c 

(vi)  monoclinic  -  1st  setting  (NSYS  =  9) 

a*  =  1/a  sin  y;  b*  =  1/b  sin  y;  c*  =  A.  a*  =  g*  =  90°;  y*  =  180°-y 

Note  y  in  the  first  setting  =  (180  -  g)  where  g  is  as  in  the  2nd  setting; 
while  a  remains  the  same  for  both  settings,  but  c  in  the  2nd  setting  becomes  b 
in  the  1st  setting  and  vice  versa. 

monoclinic  -  2nd  setting  (NSYS  »  3) 

a*  -  1/a  sin  g;  b*  =  A>  c*  =*  1/c  sin  g;  a*  »  y*  -  90°;  g  -  180°-g 

(vii)  triclinic  (NSYS  -  2) 

a*  =*  bc.sina/V;  b*  =  ca.sing/V;  c*  »  ab.sinY/V 
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8.5  Interplanar  angles  for  the  triclinic  system  in  the  direct  lattice 
between  two  planes  with  Miller  Indices  hki  and  h'k'2.1 

cos8  =  hh'a*2  +  kk'b*2  +  U'c*2  +  (hi*  +  ik' )  b*c*cosa*  n 

+  Uh'  +  hi')  c*a*cos@*  +  (hk'  +  h’k)  a*b*cot  *//(yhk4  Vk't') 

where  is  given  by 

Q.  ,  =  h2a*2  +  k2b*2  +l2c*2  +  2k£b*c*cosa*  +  2£hc*a*cos£*  +  2kha*b*cosr* 

hk4 
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9. 


APPENDIX  2 


Hie  relationship  between  the  NSYS  code3  and  Che  crystal  symmetry 


Space 

groups 

Crystal  setting  or  axes 

Point  groups 

Laue- 

symmetry 

group 

Crystal 

system 

NSYS  code  in 
FLKEFLY 

1-2 

1,  T 

T 

trlcllntc 

2 

3-15 

1st  setting 

(z-axis  unique;  GAMMA  Input) 

2,  ra,  2/ra 

2/m 

monoclinic 

9 

3-15 

2nd  setting  or  unspecified 
(y-axls  unique;  BETA  input) 

2,  ra,  2/ra 

2/m 

monocl Inlc 

J 

16-74 

222 ,  mm2  ,  mown 

minm 

orthorhombic 

4 

75-88 

4,  H  ,  4/ra 

4/m 

tetragonal 

5 

89-142 

422,  4mm,  ?2m  ,  4/mmm 

4/mmra 

tetragonal 

10 

143-148 

rhorabohedral  or  unspecified 

3,  7 

7 

trigonal 

7 

143-148 

hexagonal 

3,  7 

7 

t  rlgona 1 

6 

149-167 

rhorabohedral  or  unspecified 

32,  3m.  7m 

7m 

trigonal 

14 

149-167 

hexagonal 

32,  3m  ,  1m 

3ra 

trigonal 

1  1 

168-176 

6,  ,  6/ra 

6/m 

hexagonal 

12 

177-194 

622,  6mra,  6ra2,  6/miran 

6/ rrurun 

hexagonal 

1  3 

195-206 

23,  m3 

ra3 

cubic 

8 

207-230 

432,  53m,  m3ra 

m3no 

cubic 

15 

N0-R-9ll(S) 


21 


Page  22  blank 


10 


APPENDIX  3 


List  of  X-ray  wavelengths  pre-set  within  FIREBIRD3 
and  the  corresponding  Input/output  data  codes 


X-radiation  Corresponding  Wavelength 

WAVECO  code  applied  (A) 


Cu 

Kax 

1 

1.54050 

Cu 

Ka2 

2 

1.54434 

Cu 

Ka 

(weighted  mean) 

3 

1.54178 

Cu 

KB 

4 

1.39217 

Cr 

Kaj 

5 

2.28962 

Cr 

Ka2 

6 

2.29352 

Cr 

Ka 

(weighted  mean) 

7 

2.29092 

Fe 

Kax 

8 

1.93597 

Fe 

Ka2 

9 

1.93991 

Fe 

Ka 

(weighted  mean) 

10 

1.93728 

Co 

Kax 

11 

1.78890 

Co 

Ka2 

12 

1.79279 

Co 

Ka 

(weighted  mean) 

13 

1.79020 

Mo 

Kax 

14 

0.70926 

Mo 

Ka2 

15 

0.71354 

Ag 

Kax 

16 

0.55941 

Ag 

Ka 2 

17 

0.56381 

ND-R-9ll(S) 
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11.  APPENDIX  4 


A  full  description  of  the  operation  of  the  NSYM  codes 
NSYM  1  (which  defines  the  Bravais  lattice)  is  first  considered. 


Code  Conditions  operating 
on  h,  k  and  & 


Relevant  portion  of  short 
space  group  symbol  and  comments 


1  None 


2  h  +  k  +  l  =  2n 

3  h  +  k  =  2n 

4  h  +  l  =  2n 

5  k  +  l  =  2n 

6  h  +  k  =  2n 

k  +  l  =  2n 
h  +  l  ■  2n 

7  -h  +  k+  2.  =  3n 

8  h+k+i=3n 

9  h  -  k  +  l  =  3n 

10  h  -  k  =  3n 


P  (primitive)  or 

R  (rhombohedral  indexed  with  rhombohedral 
axes)  or 

C  (hexagonal,  as  in  1935  edition  of 

International  Tables  for  Crystal  Structure 
Determination,  Volume  I) 

I  (body  centred) 

C  (c  face-centred) 

B  (b  face-centred) 

A  (a  face-centred) 


F  (all  face  centred) 


R  (rhombohedral,  indexed  on  hexagonal 
axes  -  obverse  position) 

P  or  C  [see  above]  (hexagonal  indexed 
on  rhombohedral  axes) 

R  (rhombohedral,  indexed  on  hexagonal 
axes  -  reverse  position) 

H  (hexagonal,  triple  cell) 


ND-R-911(S) 
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Then  if  the  crystal  system  is  orthorhombic  for  0k£  (including  Okk) 
NSYM2  is  considered. 


Code  Conditions  operating  Relevant  portion  of  short 

on  h,  k  and  1  space  group  symbol  and  comments 


1  No  conditions  No  {100}  glide  planes  present 

2  k  =  2n  (P)  or  (A)  or  (I)  then  b--  or  -b- 

3  l  =  2n  (P)  or  (C)  then  c--  or  -c- 

4  k  +  l  =  2n  (P)  then  n--  or  -n- 

5  k  +  l  =  4n  (F)  then  d--  or  d- 


Then  if  the  crystal  system  is  monoclinic  (second  setting), 
orthorhombic,  or  tetragonal  or  cubic  for  h0£  (which  includes  hOh)  NSYM  3 
is  considered. 


Code 

Conditions  operating 
h,  k  and  l 

Relevant  portion  of 
space  group  symbol  and 

short 

comments 

1 

No  conditions 

No  (010)  glide  planes  present 

2 

h  =  2n 

f(P)  or  (A)  or  (I)  then  -a- 
{(P)  then  (-)-b- 

3 

l  =  2n 

(P)  or  (C)  or  (I'  then  (-)-b- 

or(-)-c-  or 

4 

h  +  l  =  2n 

(P)  then  (-)-n- 

5 

h  +  1  =  4n 

(F)  then  -d- 

Then  if  the  crystal  system  is  monoclinic  (first  setting), 
orthorhombic,  or  tetragonal  or  cubic  for  hkO  (including  hhO),  NSYF1  4  is 
considered. 


Code 

Conditions  operating 
h,  k  and  2. 

Relevant  portion  of  short 
space  group  symbol  and  comments 

1 

No  conditions 

No  (001)  glide  planes  present 

2 

h  =  2n 

(P)  or  (C)  or  (1)  --a;  or  (P--a;  or 

1  (-)a--j  or  (P)  or  (I)a 

3 

k  =  2n 

(P)  b-  or  (P)  -b 

4 

h  +  k  =  2n 

(P)--n;  or  P(-)  -n-;  or  (P)  n;  or  (P)-n-- 

5 

h  +  k  =  4n 

(F)--d;  or  (F)d 

ND-R-91 1 ( S ) 
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Then  if  the  crystal  is  tetragonal,  trigonal,  hexagonal  or  cubic  for 
hht  (which  includes  hhO  and  hhh,  but  not  OOh),  NSTM  5  is  considered. 


4 


Code 

Conditions  operating 
on  h,  k  and  8. 

Relevant  portion  of  short 
space  group  symbol  and  comments 

1 

No  conditions 

No  {110}  glide  planes  present 

2 

l  =  2n 

(P)--c;  or  (P)(-)--c;  or  (P)--n;  or  (R)  -c 

3 

2h  +  8.  =  4n 

(I)--d;  or  (I)(-)"d 

4 

h  =  2n  and  8.  =  2n 

5 

8.  =  3n 

Then 
6  is 

if  the  system  is  orthorhombic,  tetragonal  or  cubic  for  hOO,  NSYM 
considered. 

Code 

Conditions  operating 
on  h,  k  and  8. 

Relevant  portion  of  short 
space  group  symbol  and  comments 

1 

No  conditions 

No  100  screw  axes  present 

2 

h  *  2n 

(P)  2i~  or  (l)-2l;  or  (P)42 

3 

h  =  4n 

(R)  or  (F)  or  (I)  4,~;  or  (P)43~ 

Then 

OkO, 

if  the  system  is  monoclinic  (second  setting)  or  orthorhombic  for 

NSYM  7  is  considered. 

Code 

Conditions  operating 
on  h,  k  and  8. 

Relevant  portion  of  short 
space  group  symbol  and  comments 

1 

No  conditions 

No  [010]  screw  axes  present 

2 

k  =  2n 

(1)2,--;  or  ( P) ~2 3 “ 

ND-R-911(S) 
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Then  if  the  system  is  monoclinic  (first  setting),  orthorhombic, 
tetragonal,  trigonal  or  hexagonal  for  OOi.,  NSYM  8  is  considered. 


Code 

Conditions  operating 
on  h,  k  and  l 

Relevant  portion  of  short 
space  group  symbol  and  comments 

1 

No  conditions 

Neither  [001]  nor  z  screw  axes  present 

2 

l  =  2n 

(P)22  or  22-  or  42  or  42-  or  (P)--23 
(P)63  or  63-  or  63-- 

3 

l  =  3n 

( P )3 !  or  3j-b--  or  32  or  32-  or 
(P)62  or  62--  or  64  or  6„ — 

4 

l  =  4n 

(P)  or  (1)43;  or  (P)4 3 ;  or  (1)43; 
or  (P  or  I)42 — ;  or  (P)43 — 

5 

i  =  6n 

(P)62  or  6S;  or  (P)6! —  or  6S-- 

Finally  if  the  system  is 
corresponds  to  hOS.)  NSYM 

trigonal  or  hexagonal  for  hhil  (which 

9  is  considered. 

Code 

Conditions  operating 
on  h,  k  and  l 

Relevant  portion  of  short 
space  group  symbol  and  comments 

1 

No  conditions 

No  {1120}  glide  planes  present 

2 

l  =  2n 

(R)-c  or  (P)-c-  or  (P)( - )-c- 

ND-R-911(S) 
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APPENDIX  5 


A  description 

of  the 

NSYM 

and  NSYS  codes  for  each 

ot  cne 

230  space  groups 

Short 

Space  group  space 

NSYM  codes 

Shot  t 
point 

NSYS 

number  group 

symbol 

NFl 

NF2 

NF3  NF4  NFS  NF6 

NF7 

NFS  NF9  *roup 
symbol 

code 

TRICLINIC 


I 

PI 

l 

l 

l 

l 

l 

l 

l 

l 

l 

2 

2 

pT 

l 

l 

l 

l 

1 

l 

1 

1 

T 

i 

2 

MONOCLINIC 

3  (1st 

setting) 

P2 

1 

1 

l 

1 

1 

l 

l 

1 

l 

2 

9 

3  (2nd 

“  ) 

P2 

l 

1 

1 

1 

l 

1 

1 

l 

l 

2 

3 

4  (1st 

-  ) 

P21 

l 

l 

1 

l 

l 

l 

l 

2 

l 

2 

9 

4  (2nd 

•*  ) 

P21 

l 

1 

l 

1 

l 

l 

2 

l 

I 

2 

3 

5  (1st 

••  ) 

B2 

4 

1 

1 

1 

l 

1 

1 

1 

l 

2 

9 

5  (2nd 

'•  ) 

C2 

3 

1 

1 

1 

l 

l 

l 

1 

1 

2 

3 

6  (1st 

-  ) 

Pm 

1 

1 

1 

1 

l 

l 

1 

1 

l 

m 

9 

6  (2nd 

-  ) 

Pm 

1 

1 

l 

1 

1 

l 

l 

1 

1 

m 

3 

7  (1st 

) 

Pb 

1 

1 

1 

3 

1 

1 

1 

l 

l 

m 

9 

7  (2nd 

“  ) 

Pc 

1 

1 

3 

1 

l 

l 

l 

1 

1 

m 

3 

8  (1st 

-  ) 

Bm 

4 

1 

1 

l 

1 

l 

l 

1 

1 

m 

9 

8  (2nd 

•  > 

Cm 

3 

l 

1 

l 

l 

1 

1 

l 

l 

m 

3 

9  (1st 

■■  ) 

Bb 

4 

1 

1 

3 

l 

l 

1 

1 

1 

3 

9 

9  (2nd 

‘  ) 

Cc 

3 

l 

3 

l 

l 

l 

1 

l 

1 

m 

3 

10  (1st 

•’  ) 

P2/m 

1 

1 

l 

l 

1 

l 

l 

l 

1 

2/m 

9 

10  (2nd 

■  > 

P2/ra 

l 

1 

l 

1 

l 

1 

1 

l 

1 

2/m 

3 

11  (1st 

••  > 

P2,/m 

1 

1 

l 

1 

l 

1 

1 

2 

1 

2/m 

9 

ll  (2nd 

1  ) 

P2,/m 

l 

1 

1 

l 

l 

1 

2 

l 

l 

2/m 

3 

12  (1st 

••  > 

B2/m 

4 

l 

l 

l 

l 

l 

l 

1 

1 

2/m 

9 

12  (2nd 

•’  ) 

C2/m 

3 

1 

l 

1 

l 

l 

1 

l 

1 

2/m 

3 

13  (Isc 

•*  ) 

P2/b 

1 

l 

l 

3 

l 

1 

l 

l 

l 

2/m 

9 

13  (2nd 

■*  ) 

P2/c 

l 

l 

3 

1 

1 

1 

1 

l 

1 

2/m 

3 

14  (1st 

■  > 

P2,/b 

1 

l 

l 

3 

l 

1 

1 

2 

1 

2/m 

9 

14  (2nd 

•  ) 

P2,/c 

1 

l 

3 

l 

l 

1 

2 

1 

1 

2/m 

3 

15  (Isc 

•*  ) 

827b 

4 

1 

1 

3 

l 

1 

1 

1 

1 

2/m 

9 

15  (2nd 

> 

C2/c 

3 

l 

3 

1 

1 

l 

l 

l 

1 

2/m 

3 

ORTHORHOMBIC 

16 

P222 

l 

l 

l 

l 

1 

1 

l 

l 

l 

222 

4 

17 

P2221 

l 

l 

1 

l 

l 

l 

l 

2 

1 

222 

4 

18 

P2  2  2 

1 

l 

1 

l 

1 

2 

2 

1 

1 

222 

4 

19 

P2,272, 

1 

1 

1 

1 

l 

2 

2 

2 

1 

222 

4 

20 

C2$2l 

3 

l 

1 

l 

l 

l 

l 

2 

1 

222 

4 

21 

C222 

3 

l 

l 

l 

l 

1 

1 

1 

I 

222 

4 

22 

F222 

6 

l 

1 

l 

1 

l 

l 

1 

1 

222 

4 

23 

1222 

2 

l 

l 

l 

l 

l 

1 

1 

l 

222 

4 

24 

12,2,2, 

2 

l 

1 

1 

l 

1 

1 

1 

l 

222 

4 

25 

Pran2 

1 

l 

l 

l 

1 

l 

1 

l 

1 

mm2 

4 

26 

Pmc2L 

l 

1 

3 

1 

1 

1 

l 

1 

l 

mm2 

4 

27 

Pcc2 

l 

3 

3 

l 

l 

t 

l 

1 

l 

mm? 

4 

28 

Pma2 

1 

l 

2 

1 

l 

l 

1 

l 

l 

mm2 

4 

29 

Pca21 

1 

3 

2 

1 

I 

L 

l 

1 

l 

mm2 

4 

30 

Pnc2 

1 

4 

3 

1 

l 

l 

1 

l 

l 

mm2 

4 

31 

Pmn2. 

l 

1 

4 

L 

1 

1 

l 

1 

l 

mm? 

4 

32 

Pba2 

l 

2 

2 

1 

l 

l 

l 

l 

1 

mm2 

4 

33 

Pna2 

1 

4 

2 

1 

l 

1 

l 

2 

l 

mm2 

4 

34 

Pnn2 

l 

4 

4 

1 

1 

1 

l 

1 

1 

mm2 

4 

35 

Can2 

3 

l 

1 

l 

l 

l 

l 

l 

l 

mm2 

4 

36 

C*c2, 

Ccc2* 

3 

l 

3 

1 

l 

1 

l 

1 

1 

mm2 

4 

37 

3 

3 

3 

l 

1 

1 

l 

1 

1 

mm2 

4 

38 

Aan2 

5 

l 

1 

1 

l 

1 

1 

l 

1 

mm2 

4 

39 

Abn2 

5 

2 

l 

1 

l 

1 

1 

» 

1 

mm2 

4 

40 

Ama2 

5 

l 

2 

l 

l 

l 

1 

1 

l 

mm2 

4 

41 

Aba2 

5 

2 

2 

1 

L 

1 

1 

1 

l 

mm2 

42 

Fm»2 

6 

l 

1 

l 

l 

l 

L 

1 

L 

mm2 

4 

43 

Fdd2 

6 

5 

5 

1 

1 

1 

l 

l 

1 

mm2 

4 

44 

Imm2 

2 

l 

l 

l 

i 

l 

l 

l 

l 

mm2 

4 

45 

Iba2 

2 

2 

2 

1 

1 

1 

l 

l 

1 

™2 

4 

MD-R-9LKS) 
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Continued 


Appendix  5  continued 


Space  group 

Short 

space 

NSW  codes 

Short 

point 

NSYS 

number 

group 

symbol 

NFl 

NF2 

NF3 

NF4  NFS  NF6 

NF7 

NF8 

NF9 

group 

symbol 

code 

ORTHORHOMBIC 

46 

(continued) 

Iaa2 

2 

1 

2 

1 

1 

l 

1 

1 

L  mm? 

4 

47 

Pmom 

l 

1 

1 

1 

1 

l 

l 

1  : 

l  ama 

4 

48 

Pnnn 

l 

4 

4 

4 

1 

1 

1 

1 

l  BOMB 

4 

49 

Pccm 

l 

3 

3 

l 

l 

l 

l 

1  ] 

mas 

4 

50 

Pban 

1 

2 

2 

4 

l 

1 

l 

1  ! 

L  ma. 

4 

51 

?mu 

l 

l 

l 

2 

1 

1 

l 

1  1 

_ 

4 

52 

Pnna 

l 

4 

4 

2 

1 

1 

l 

1  1 

mas 

4 

53 

Pmna 

1 

1 

4 

2 

l 

1 

l 

1  1 

mmm 

4 

54 

Pcca 

1 

3 

3 

2 

1 

1 

l 

l  1 

mo a 

4 

55 

Pbam 

1 

2 

2 

l 

l 

l 

l 

l  1 

Him 

4 

56 

Pccn 

1 

3 

3 

4 

1 

l 

l 

l  1 

mian 

4 

57 

Pbca 

1 

2 

3 

l 

1 

1 

1 

1  ] 

Timm 

4 

58 

Pnna 

1 

4 

4. 

1 

1 

l 

l 

1  1 

mmm 

4 

59 

Prnum 

1 

1 

1 

4 

1 

1 

1 

1  1 

mum 

4 

60 

Pbcn 

1 

2 

3 

4 

1 

1 

1 

1  J 

!  » 

4 

61 

Pbca 

l 

2 

3 

2 

l 

1 

1 

1  1 

man 

4 

62 

Pnma 

1 

4 

l 

2 

l 

1 

l 

t  1 

mum 

4 

83 

Coca 

3 

1 

3 

1 

1 

l 

1 

1  I 

mmm 

4 

64 

Coca 

3 

l 

3 

2 

l 

l 

l 

1  I 

man 

4 

65 

Cflfflfli 

3 

1 

l 

1 

1 

1 

1 

I  1 

oam 

4 

66 

Cccm 

3 

3 

3 

1 

l 

l 

l 

1  1 

nraa 

4 

67 

Came 

3 

1 

l 

2 

1 

l 

l 

l  1 

man 

4 

68 

Ccca 

3 

3 

3 

2 

1 

1 

l 

l  1 

mmm 

4 

69 

Fmiaa 

6 

1 

l 

1 

1 

1 

1 

1  1 

mrani 

4 

70 

Fddd 

6 

5 

5 

5 

1 

1 

1 

1  1 

mmm 

4 

71 

Tmrm 

2 

l 

l 

l 

1 

1 

1 

1  : 

, 

4 

72 

Ibara 

2 

2 

2 

l 

l 

l 

1 

1 

L  asms 

4 

73 

Ibca 

2 

2 

3 

2 

l 

l 

1 

l  ! 

1  mnva 

4 

74 

bnu 

2 

1 

l 

2 

l 

1 

1 

1 

,  mmm 

4 

TETRAGONAL 

75 

P4 

1 

1 

1 

1 

l 

1 

l 

t 

l  4 

5 

76 

P4, 

l 

1 

1 

L 

l 

l 

1 

4 

l  4 

5 

77 

«2 

l 

1 

l 

l 

1 

l 

l 

2 

L  4 

5 

78 

P43 

l 

l 

l 

l 

1 

l 

l 

4 

l  4 

5 

79 

14 

2 

1 

l 

1 

l 

l 

l 

1 

l  4 

5 

80 

14, 

2 

l 

l 

l 

l 

l 

l 

4  : 

L  4 

5 

81 

P4l 

1 

J 

l 

1 

1 

1 

l 

1 

l  4 

5 

82 

14 

2 

1 

l 

1 

1 

l 

l 

1  : 

I  4 

5 

33 

P4/a 

1 

1 

1 

1 

l 

l 

1 

1 

1  4/m 

5 

84 

P42/m 

1 

1 

l 

1 

1 

1 

1 

2 

1  4/m 

5 

85 

P4/n 

1 

l 

l 

4 

1 

l 

l 

1 

1  4/o 

5 

86 

P42/n 

I4/a 

1 

I 

1 

4 

1 

1 

1 

2 

1  4/m 

5 

87 

2 

l 

t 

1 

1 

l 

l 

1 

l  4/m 

5 

88 

14 1 /  a 
P422 

2 

1 

l 

2 

l 

l 

1 

4 

1  4/m 

5 

89 

1 

l 

l 

l 

1 

1 

l 

1 

l  422 

to 

90 

P42.2 

l 

l 

1 

1 

1 

2 

1 

1 

l  422 

10 

91 

P4.i2 

1 

1 

l 

l 

1 

1 

1 

4 

l  422 

10 

92 

P4j2  2 

l 

1 

l 

l 

l 

2 

l 

4 

l  422 

10 

93 

1 

1 

1 

1 

1 

1 

1 

2 

1  422 

10 

94 

P4z2l2 

1 

1 

1 

1 

l 

2 

l 

2 

l  422 

10 

95 

P4322 

1 

1 

l 

1 

l 

1 

1 

4 

1  422 

10 

96 

P4.2.2 

l 

1 

l 

1 

1 

2 

1 

a 

l  422 

10 

97 

14I21 

2 

l 

1 

l 

l 

1 

1 

l 

l  422 

10 

98 

I4t22 

2 

l 

l 

t 

1 

1 

l 

4 

l  422 

10 

99 

P4oa 

1 

l 

l 

1 

l 

l 

l 

1 

1  4mm 

10 

100 

P4ba 

1 

1 

2 

l 

1 

l 

l 

l 

1  4(nB 

10 

101 

P42cm 

l 

l 

3 

l 

l 

1 

L 

L 

[  4mm 

10 

102 

P42na 

P4cc 

1 

l 

4 

l 

l 

l 

1 

1 

1  4  mrn 

10 

103 

l 

1 

3 

l 

2 

1 

1 

l 

1  4mm 

10 

104 

P4nc 

l 

1 

4 

1 

2 

l 

1 

l 

1  inm 

10 
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Appendix  5  continued 


Short 

NSW 

codes 

Short 

Space  group 

space 

point 

NSYS 

number 

group 

NF1 

NF2 

NF3 

NF4 

NFS 

NF6 

NF7 

NF8 

NF9 

group 

code 

synboi 

symbol 

TETRAGONAL  (Continued) 

105 

P42ac 

1 

l 

1 

l 

2 

l 

l 

I 

1 

4oo 

10 

LQ6 

P42bc 

1 

1 

2 

l 

2 

1 

l 

1 

1 

4  an 

10 

107 

I4oa 

2 

l 

1 

1 

l 

l 

l 

1 

1 

4  an 

10 

108 

I4cm 

2 

1 

3 

1 

l 

l 

l 

1 

1 

4  an 

10 

109 

14  ^md 

2 

1 

l 

1 

3 

l 

l 

1 

1 

4mn 

10 

110 

111 

14  .cd 
?Zim 

2 

1 

l 

1 

3 

l 

1 

l 

3 

1 

1 

l 

1 

1 

1 

1 

1 

1 

4«m 

42m 

10 

10 

112 

P42c 

1 

1 

1 

l 

2 

l 

1 

1 

1 

42m 

10 

113 

P42 .  m 

l 

1 

1 

1 

l 

2 

L 

1 

l 

42m 

10 

114 

?42rc 

1 

l 

1 

l 

2 

2 

l 

I 

1 

42m 

10 

115 

P7m2 

1 

l 

l 

1 

l 

l 

l 

1 

l 

42m 

10 

116 

P4c2 

1 

1 

3 

l 

1 

1 

1 

1 

l 

42m 

10 

117 

Plt>2 

1 

1 

2 

l 

1 

l 

1 

l 

1 

42m 

L0 

118 

P4n2 

l 

1 

4 

l 

1 

l 

l 

1 

1 

42m 

10 

119 

I4o2 

2 

1 

1 

l 

l 

l 

l 

1 

l 

_42m 

10 

120 

I4c2 

2 

1 

3 

1 

1 

1 

l 

1 

1 

T2m 

10 

121 

I42m 

2 

1 

1 

l 

1 

1 

1 

1 

1 

42m 

10 

122 

I42d 

2 

l 

l 

1 

3 

1 

l 

1 

1 

42a 

10 

123 

P/4nanm 

1 

l 

l 

1 

1 

l 

L 

L 

l 

4 1  tsmm 

10 

124 

P4/mcc 

l 

l 

3 

t 

2 

l 

l 

1 

1 

4/tnnw 

10 

125 

P4/ nbra 

l 

2 

1 

4 

1 

l 

l 

1 

l 

4  /  nunra 

10 

126 

P4/nnc 

l 

4 

1 

4 

2 

l 

l 

1 

l 

4  /  mm 

10 

127 

P4/obm 

l 

2 

1 

l 

1 

2 

l 

1 

1 

4/mrmn 

10 

128 

P4/nmc 

1 

4 

1 

l 

2 

2 

l 

1 

l 

4  /  nunm 

10 

129 

P4/  mam 

1 

1 

1 

4 

l 

2 

1 

1 

L 

4  /  Hiram 

10 

130 

P  4/tvcc 

l 

3 

1 

4 

2 

2 

1 

2 

1 

4/rmmn 

10 

131 

P42/osac 

1 

1 

l 

1 

2 

1 

1 

2 

1 

4  /  mmm 

10 

132 

P42/mca 

1 

3 

1 

1 

1 

1 

l 

2 

1 

4 /aura 

10 

133 

P42/nbc 

1 

2 

1 

4 

2 

2 

1 

2 

L 

4  /  nunm 

10 

134 

P42/nnm 

1 

4 

1 

4 

1 

2 

l 

2 

1 

4/omra 

10 

135 

p42/(abc 

l 

2 

1 

l 

2 

2 

l 

2 

1 

4  / Hiram 

10 

136 

P42/mnm 

l 

4 

l 

l 

1 

2 

l 

2 

l 

4  /  mmm 

10 

137 

p42/nac 

l 

i 

1 

4 

2 

2 

1 

2 

l 

4/aanm 

10 

138 

P42/ncm 

14 /anno 

1 

3 

1 

4 

l 

2 

1 

2 

1 

4  /mmm 

10 

139 

2 

1 

i 

1 

l 

2 

l 

2 

l 

4  /  mmm 

10 

140 

14  /mem 

2 

3 

1 

l 

1 

2 

1 

2 

l 

4/nmm 

10 

141  #* 

14. /aod 

2 

3 

1 

2 

3 

2 

l 

4 

1 

4/tnmm 

10 

142*# 

14^/acd 

2 

3 

X 

2 

3 

2 

1 

4 

1 

4/nmm 

10 

TRIGONAL 

143 

P3 

1 

l 

1 

l 

l 

1 

l 

1 

l 

3 

b 

144 

1 

1 

1 

l 

1 

1 

l 

3 

1 

3 

6 

145 

P32 

1 

1 

1 

l 

l 

l 

l 

3 

l 

3 

6 

146  (rhombohedrai  axes) 

R32 

l 

l 

l 

l 

1 

i 

1 

1 

1 

3 

7 

146  (hexagonal  axes) 

R3 

7 

1 

1 

I 

5 

L 

1 

l 

l 

3 

6 

147 

P T 

l 

1 

1 

1 

l 

1 

1 

1 

1 

6 

148  (rhombohedrai  axes) 

R3 

l 

l 

1 

1 

l 

1 

l 

l 

l 

3 

7 

148  (hexagonal  axes) 

R3 

7 

l 

1 

l 

5 

1 

1 

1 

l 

3 

6 

149 

P312 

1 

1 

1 

l 

l 

l 

l 

l 

1 

312 

11 

150 

P321 

1 

1 

1 

1 

l 

1 

l 

1 

l 

321 

11 

151 

P3i  12 

l 

1 

1 

l 

l 

l 

l 

3 

l 

312 

11 

152 

P3l21 

1 

l 

•  l 

1 

1 

l 

l 

3 

l 

321 

11 

153 

P3212 

1 

l 

1 

1 

l 

1 

l 

3 

l 

312 

LI 

154 

P3221 

R32 

1 

l 

1 

l 

l 

l 

l 

3 

1 

321 

11 

155  (rhombohedrai  axes) 

l 

1 

l 

l 

l 

1 

l 

1 

l 

32 

14 

155  (hexagonai  axes) 

R32 

7 

1 

l 

1 

5 

l 

l 

1 

1 

32 

11 

156 

P30! 

1 

l 

l 

1 

l 

1 

1 

1 

i 

3m  1 

l  1 

157 

P3  Ira 

l 

l 

l 

1 

1 

l 

1 

1 

l 

31m 

ll 

158 

P3cl 

l 

1 

L 

l 

l 

l 

1 

2 

2 

3  m  1 

11 

159 

P31c 

1 

l 

1 

l 

2 

l 

l 

2 

L 

31m 

Ll 
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Appendix  5  continued 


Space  group 
number 

Short 

space 

group 

symbol 

NF1 

NF2 

NSYM 

NF3  NF4 

codes 

NFS 

NFo 

NF7 

NF8 

NF9 

Short 
point 
group 
syraoo  1 

NSYS 

code 

TRIGONAL  (continued) 

160  ( rhombohedral  axes) 

R3m 

1 

1 

l 

1 

1 

1 

1 

1 

1 

3m 

1- 

160  (hexagonal  axes) 

R3m 

7 

l 

1 

1 

5 

1 

l 

l 

1 

3m 

1  1 

161  (rhombohedral  axes) 

R3c 

1 

I 

1 

l 

4 

1 

1 

1 

l 

3m 

1  - 

1 6 L  (hexagonal  axes) 

R3c 

7 

1 

1 

l 

5 

1 

1 

5 

2 

3m 

1 1 

162 

rtu 

1 

1 

1 

l 

l 

1 

1 

1 

1 

31m 

1  1 

163 

P31c 

l 

1 

l 

1 

2 

l 

l 

2 

1 

31m 

1  1 

164 

P3ml 

1 

l 

1 

l 

1 

1 

1 

1 

1 

3mi 

1  1 

165 

P3cl 

l 

1 

l 

1 

l 

1 

1 

z 

2 

3  ml 

L  1 

166  (rhoobohedrai  axes) 

R3m 

1 

l 

1 

1 

1 

1 

1 

1 

1 

3m 

14 

166  (hexagonal  axes) 

R3m 

7 

l 

1 

1 

5 

1 

1 

1 

1 

3m 

1  1 

167  (rhombohedral  axes) 

Rlc 

1 

1 

l 

1 

2 

1 

L 

1 

1 

To 

167  (hexagonal  axes) 

R3c 

7 

1 

1 

1 

5 

1 

1 

5 

2 

3m 

I  i 

HEXAGONAL 

168 

P6 

1 

1 

l 

1 

1 

1 

1 

1 

1 

■ 

169 

P6i 

1 

1 

l 

1 

1 

1 

1 

5 

1 

0 

i_2 

170 

P65 

1 

1 

1 

1 

1 

1 

1 

5 

1 

o 

12 

171 

P69 

1 

1 

1 

1 

1 

1 

1 

3 

l 

6 

14 

172 

1 

1 

1 

1 

1 

1 

1 

3 

1 

0 

1- 

173 

l 

1 

l 

1 

1 

1 

1 

2 

1 

6 

iJ 

174 

P6 

1 

1 

1 

l 

1 

1 

1 

1 

1 

0 

L2 

175 

P6/ra 

1 

1 

l 

1 

1 

l 

1 

1 

1 

6/  m 

12 

176 

P63/m 

P622 

1 

1 

1 

1 

1 

l 

l 

2 

1 

o/m 

1. 

177 

1 

1 

1 

l 

1 

1 

l 

1 

1 

622 

:  3 

178 

P6.22 

1 

1 

1 

1 

1 

i 

1 

5 

1 

o22 

13 

179 

P6c22 

P6222 

l 

1 

1 

1 

1 

1 

l 

5 

1 

622 

1  > 

180 

1 

1 

l 

l 

1 

1 

1 

3 

1 

oZ2 

i  i 

181 

P6422 

1 

1 

1 

1 

1 

1 

1 

3 

l 

o22 

l  3 

182 

P6322 

1 

1 

l 

1 

1 

1 

1 

2 

1 

622 

1  >■ 

183 

P6mm 

1 

l 

l 

1 

1 

l 

l 

1 

1 

omo 

.  1 

184 

P6cc 

1 

l 

l 

l 

2 

l 

l 

2 

2 

oram 

i  3 

185 

Pb^cm 

P6_3mc 

1 

1 

1 

1 

l 

1 

1 

2 

2 

omm 

1  3 

186 

1 

1 

l 

1 

2 

l 

l 

2 

1 

omm 

1  3 

187 

P6ra2 

1 

1 

1 

1 

1 

1 

l 

1 

1 

oraJ 

1  3 

188 

P5c2 

1 

1 

1 

l 

1 

1 

1 

2 

2 

"braJ 

13 

189 

P62ra 

1 

1 

1 

l 

1 

1 

1 

1 

1 

&2m 

i  i 

190 

P62c 

1 

l 

1 

1 

2 

1 

1 

2 

1 

62m 

i  3 

191 

P6/naan 

1 

1 

1 

l 

l 

1 

1 

1 

l 

o  /  romro 

1  3 

192 

P6/mcc 

1 

1 

l 

l 

2 

l 

I 

2 

2 

r*/  mmm 

13 

193 

Pb^/mcm 

P63/camc 

1 

l 

1 

1 

1 

l 

1 

2 

2 

o  /  mmE 

i  i 

194 

1 

1 

1 

l 

2 

1 

1 

2 

l 

inmra 

I  3 

CUBIC 

195 

P23 

1 

1 

1 

1 

l 

t 

1 

l 

1 

2  3 

b 

196 

F23 

6 

1 

1 

1 

l 

l 

1 

1 

1 

23 

0 

197 

123 

2 

1 

1 

1 

1 

1 

l 

1 

1 

23 

» 

198 

P2.3 

1 

1 

l 

1 

1 

2 

1 

J 

l 

23 

0 

199 

12,3 

2 

1 

1 

1 

1 

1 

1 

1 

1 

23 

0 

200 

Pm3 

l 

l 

1 

1 

1 

l 

1 

1 

1 

m3 

» 

201 

Pn3 

1 

l 

1 

4 

1 

2 

1 

1 

1 

m3 

8 

202 

Foil 

6 

1 

l 

1 

l 

1 

1 

1 

1 

m3 

a 

203 

Fdl 

6 

l 

1 

5 

1 

3 

1 

1 

1 

m3 

0 

204 

lm3 

2 

1 

1 

1 

l 

l 

1 

1 

1 

m3 

b 
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Appendix  3  continued 


Space  group 

Shore 

space 

NSYM  codes 

Short 

point 

l.S'tS 

number 

group 

symbol 

NF1 

NF2 

HF3 

NF4  NFS  NFb 

NF7  NFS 

NF9 

group 

symbol 

code 

CUBIC  (continued) 

205 

Pa3 

l  l 

1  2 

l 

2  l 

l  l  m3 

tt 

20b 

U3 

2  i 

l  2 

1 

l  l 

1  1  m3 

6 

207 

P432 

1  l 

l  l 

1 

l  1 

l  l  432 

15 

206 

P4,32 

l  1 

l  1 

i 

2  l  l 

l  1  432 

15 

209 

F4l2 

6  1 

l  1 

1 

1  I 

l  l  432 

15 

210 

F4.32 

6  1 

l  L 

I 

3  l 

l  I  432 

15 

211 

1432 

2  l 

l  l 

1 

1  1 

l  1  432 

15 

212 

P4332 

P4.32 

l  L 

I  l 

1 

3  1 

l  L  432 

15 

213 

l  l 

l  1 

l 

3  t 

l  L  432 

15 

214 

14  v  32 

2  l 

l  l 

l 

3  l 

L  l  432 

15 

213 

P43n 

I  l 

l  t 

l 

l  l 

L  t  -3ra 

13 

216 

F4  3ra 

6  l 

L  I 

l 

1  1 

L  1  -3m 

13 

217 

143  m 

2  I 

1  L 

1 

1  I 

l  1  -3m 

15 

218 

P43n 

l  l 

l  l 

2 

2  I 

l  1  4  3m 

15 

219 

F43c 

6  l 

l  1 

2 

1  I 

l  1  -3ra 

15 

220 

14  3d 

2  I 

1  l 

3 

3  i 

l  l  -'3m 

15 

221 

PttJn 

1  1 

l  l 

1 

1  l 

l  l  mini 

15 

222 

Pn3 

l  l 

1  4 

2 

2  1 

L  1  mlm 

15 

223 

Pm3n 

l  l 

1  1 

2 

2  l 

l  l  m3m 

15 

224 

Pn3ra 

l  l 

1  1 

l 

2  l 

l  l  mlm 

15 

225 

Fm3m 

6  l 

I  \ 

l 

l  l 

1  l  mini 

15 

226 

Ffuic 

6  1 

l  1 

2 

I  l 

l  l  mlm 

15 

227 

Fd3ra 

6  1 

1  5 

l 

3  1 

l  I  mini 

1  5 

228 

Fd3c 

6  l 

l  5 

2 

3  l 

L  l  ralrn 

15 

229 

Ira3m 

2  I 

l  1 

1 

l  l 

l  l  mlm 

l  5 

230 

laid 

2  1 

1  2 

3 

3  l 

l  1  '  ralm 

15 
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13.  APPENDIX  6 


Role  of  the  various  subroutines  used  in  FIREBIRD3 


File  No. 


Operation 


S/R  FBI  Main  control  and  starting  routine. 

S/R  FB2  To  read  in  data  off  the  parameter  cards. 

S/R  FB3  To  calculate  the  values  of  the  a*  b*  c*  a*  etc. 

S/R  FB4  To  output  and  calculate  the  parameters  and  dimensions  of  the 


S/R  FB5  To  write  the  headings  for  the.  output  of  the  Bragg  angles. 

S/R  FB6  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 

reflections  for  the  cube. 

S/R  FB7  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 
reflections  for  the  triclinic  case. 

S/R  FB8  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 
reflections  for  the  trigonal/rhombohedral  case. 

S/R  FB9  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 
reflections  for  the  trigonal/hexagonal  case. 

S/R  FB10  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 
reflections  for  the  monoclinic  case. 

S/R  FB11  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 
reflections  for  the  orthorhombic  case. 

S/R  FB12  To  work  out  the  Bragg  angles  for  the  hexagonal  crystal  system. 

Relevant  HKL  combinations  are  created  and  their  validity  is 
then  checked,  if  required,  by  S/R  FB19. 

S/R  FB13  To  work  out  the  Bragg  angles  and  to  eliminate  unnecessary 
reflections  for  the  tetragonal  case. 

S/R  FB14  To  sort  out  the  Bragg  angles  into  ascending  order. 

S/R  FB15  To  output  the  calculated  and  observed  Bragg  angles. 

S/R  FB16  To  output  a  table  of  hkl's  for  input  to  the  program  "FIREFLY". 

S/R  FB17  To  compare  the  reference  plane  Miller  indice-  and  to  calculate 

the  interplanar  angle. 

S/R  FB18  To  modify  and  produce  combinations  of  unit  cell  parameter 
values  for  each  crystal  system.  Parameters  that  may  be 
modified  depend  on  the  crystal  system,  but  are  modified 
independently  of  each  other  by  user  supplied  increments  until 
the  limiting  value  (user  supplied)  is  reached. 
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S/R  FB19  To  eliminate  unwanted  hkl's  using  given  conditions,  defined  by 
the  NSYM  codes. 

S/R  FB20  To  identify  the  wavelengths  using  the  given  wave  code. 

S/R  FB21  To  determine  the  dimensions  of  the  crystal,  which  have  not  been 

supplied. 

S/R  FB22  To  determine  the  value  of  NSYS  for  given  values  of  SG  and 
SETT. 

S/R  FB23  Function  used  in  the  S/R  FB24. 

S/R  FB24  Infers  the  multiplicity  of  a  reflection  on  the  basis  of  the 
Laue  group  and  its  Miller  indices. 

S/R  FB25  To  create  hkl  values  of  ascending  order,  k  may  not  be  >  h,  and 
1  may  not  be  >  k. 

S/R  FB26  To  calculate  the  limiting  h  value  for  a  crystal  system. 

3/R  FB27  Retrieves  results  which  have  been  output  to  disc  by  S/R  FB15. 
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14.  APPENDIX  7 

Instructions  for  compiling  and  running  FIREBIRD  3  on 
the  UKAEA's  Rislev  ICL  3980  Computer 

1 .  The  Program  Modules 

FIREBIRD  comprises  a  main  program,  26  subroutines  and  1  function. 

For  convenience  these  routines  are  stored  in  7  files  in  the  library 
: RFLBOF . FIREBIRDSCE .  These  are: 

FIRE1,  FIRE2,  FIRE3,  FIRE4,  FIRE5,  FIRE6  and  FIRE 7 

The  main  program  and  subroutines  FB2  and  FB3  are  stored  in  FIRE1 . 
Subroutines  FB4,  FB5  and  FB6  are  stored  in  FIRE2. 

Subroutines  FB7,  FB8  and  FB9  are  stored  in  FIRE3. 

Subroutines  FB10,  FB11,  FB12  and  FB13  are  stored  in  FIRE4. 

Subroutines  FB14,  FB15,  FB16  and  FB17  are  stored  in  FIRE5. 

Subroutines  FB18,  FB19  and  FB20  are  stored  in  FIRE6. 

Subroutines  FB21,  FB22,  FB23,  FB24,  FB25,  FB26  and  FB27  and  the  function 
ASIN  are  stored  in  FIRE7. 

2.  Compiling  the  program 

If  any  of  the  program  subroutines  are  modified,  then  the  file(s) 
containing  the  modified  subroutines  will  have  to  be  compiled. 

For  example,  suppose  FB5  and  FB15  are  modified,  then  the  following 
commands,  typed  in  during  a  MAC  session  of  user  :RFLB0F,  will  compile  the 
program. 

PRO ( F IREBIRD ) 

C0PT( NONE , MESSAGES=NOCOMMENTS ) 

FTN(C,FIRE2) 

FTN(C,FIRE5) 

COLLECT ( COLLSCL ) 

Any  compiler  comments  at  the  end  of  the  compilation  can  be  ignored. 

For  details  on  the  ICL  3980  commands  used,  refer  to  the  3980  Users'  Manual. 

3.  Running  the  program 

The  following  file  is  needed  to  run  FIREBIRD  on  the  ICL  3980, 

JOB( :RFLB0F. jobname,ACC=account  number) 

DATA 

W0R( 9, L, 100, 500) 

ASS ( : RFLBOF . FLINTDATA ,28) 

OUT(MAX=*3000) 

R0PT( TRIES- 100 , CONT-YES ) 

USING( COLLOMF ) 

ENTER  FIREBIRD 

data 

ENDJOB 

**** 


ND-R-911(S) 


37 


If  the  data  are  to  be  read  from  a  data  file,  e.g.  the  test  data  file 
FIREBIRDSCE.FBDAT,  then  the  DATA  card  must  be  replaced  by  the  following 
card. 


DATA(FILE=FIREBIRDSCE. FBDAT) 

No  data  must  be  input  after  the  ENTER  command  if  the  above  card  is  used. 
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15.  APPENDIX  8 


Programs  in  SL's  HELLFIRE 

suite  of  programs 

(intended  for  use  in  connection  with 

X-ray  powder  dif fractometrv 

HELLFIRE1 

EPIMETHEUS( 14) 

Validates  experimental  data 

HELLFIRE2 

HELLFIRE2 

HELLFIRE2A 

FIRECRACKER(13) 
FIRECRACKER2( 15) 
FIRECRACKER3 ( 22 ) 

Calibration 
of  X-ray 
diffractometers 

HELLFIRE3A 

HELLFIRE4 

FIRESTAR2(23) 

FIRERICE  &  FIRECHOP(24) 

Determination  of  lattice 
parameters. 

Rachinger  separation  of 
ctl  and  a2  Bragg  reflections 

HELLFIRE5 

FXREBRAND(25) 

Locates  Bragg  maximum 

HELL FIRE 6 

FIREPL0T(26) 

Plots  out  diffraction 
patterns  and  can  smooth. 

HELLFIRE7 

FIREBIRD3(21) 

Calculates  the  positions  of 
Bragg  reflections. 

HELLFIRE8 

HELLFIRE8A 

HELLFIRE8B 

FIREFLY3A  (FIREBEE)(1) 
FIREFLY4( 1 ) 

FIREFLY6  (FIRECOMET)20 

Calculations  of  the 
intensities  of  Bragg 
reflections  given  the  unit 
cell  and  if  required 
the  necessary  Miller 
indices. 

HELLFIRE9 

FIREWOLF(27 ) 

Integrates  Bragg  reflection 
intensity. 

HELLFIRE10 

FENRIS(28) 

Derives  crystallite  sizes 
and  strains  from  integral 
breadths . 

HELLFIRE11 

FAFNIR(29) 

Derives  crystallite  sizes 
and  strains  from  integral 
breadths,  but  includes 
possibility  to  calculate 
Scherrer  constants  for 
various  types  of  prism. 

HELLFIRE12 

FIRENET(30) 

Plots  direct  pole  figures. 

Page  40  blank 
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